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We study the renormalized energy-momentum tensor (EMT) of cosmological scalar fluctuations 
during the slow-rollover regime for chaotic inflation with a quadratic potential and find that it is 
characterized by a negative energy density which grows during slow-rollover. We also approach the 
back-reaction problem as a second-order calculation in perturbation theory finding no evidence that 
■ the back- reaction of cosmological fluctuations is a gauge artifact. In agreement with the results for 

the EMT, the average expansion rate is decreased by the back-reaction of cosmological fluctuations. 



o 
o 



INTRODUCTION 



CO ■ After the latest CMB data , inflation (see 0, for a textbook review) seems the most promising theory in 
explaining the large scale structure of the Universe. According to inflation, the large scale structure of the Universe 
exhibits the fingerprint of quantum fluctuations amplified during the accelerated era ■ On modelling inflation and 
its transition to the standard Big Bang cosmology the constraints on the amplitude and spectrum of CMB fluctuation 
become constraints on the physics of inflation. All this modelling depends on the linear treatment of cosmological 
■ perturbations. Within the framework of inflationary cosmology, on using the most recent data, we are therefore close 
CO | to measuring the details of the spectrum of fluctuations, while not having a definite idea on the energy content of 
these fluctuations and how they back-react on the inflationary expansion responsible for their amplification. 
| Although cosmological linear perturbation theory during inflation is almost a textbook subject, the understanding 
q-i of nonlinear effects is still at the forefront of research. The nonlinearity and gauge invariance of General Relativity are 
tenacious obstacles both at technical and interpretational levels, nonetheless the gauge issue has been solved to higher 
orders in perturbation theory 0,0- Further, many interesting effects are appreciable only beyond the linear order. 
From the theoretical point of view the back-reaction of gravitational fluctuations on the geometry is one of the most 
.£h ! interesting issues Q- Within the inflationary context, this problem has been tackled by Abramo, Brandenberger and 
Mukhanov 0J • The intringuing result that the energy- momentum tensor (henceforth EMT) of fluctuations may 
slow down inflation 0|| has subsequently generated renewed interest in the subject of back-reaction 0,0,^2- The 
final answer on the physical significance of this back-reaction effect is still under debate [ToL ITU . The problem of 
gravitational back-reaction for black-holes has also been tackled for gravitational waves [l3| . More recently these non- 
linear effects have also drawn attention in connection with observations, since non-linear cosmological perturbations 
introduce non gaussian signatures in the power spectrum 0, . 

The aim of this paper is to compute the renormalized EMT of cosmological fluctuations during inflation, according 
to the adiabatic regularization scheme also used in our previous paper [17]. The model considered here is the 
slow-rollover regime of inflation driven by a massive inflaton, but we believe that the results obtained here also hold 
for other inflationary models. We find that the averaged (with respect to the adiabatic vacuum) renormalized EMT 
of cosmological fluctuations during slow-rollover is characterized by a negative energy density and a de Sitter-like 
equation of state (this result was found for long- wavelength modes in In a naive approach this would lead 
us to think that the EMT of cosmological fluctuations slows down inflation. We also evaluate the back-reaction on 
the geometry in a systematic way by proceeding self-consistently to second-order in perturbation theory. In order 
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to do this we give a systematic treatment of second-order perturbation theory for single scalar field driven inflation 
(see also 0, f° r the second-order formalism) . The gauge used for scalar perturbations in this paper is the 

uniform curvature gauge j2£j (UCG) generalized to second order. In the UCG the spatial sections are not perturbed 
by scalar fluctuations. We believe that this gauge is more convenient for our problem than the more frequently used 
longitudinal gauge. 

In this paper we focus on the EMT of scalar cosmological perturbations. This is indeed the relevant effect since 
vector perturbations decay and gravitational waves are described by an EMT which is equivalent to a massless field 
|2*lj . whose main effect is non-leading with respect to scalar fluctuations (see however Q for the two- loop calculation). 

The plan of the paper is as follow. In section 2 we present the linear cosmological perturbations in the UCG. In 
section 3 we extend the UCG to second order and we give expressions for the Einstein and the energy-momentum 
tensors. In section 4 we present and illustrate the use of the Einstein equations to second order. In sections 5 and 7 
we give the approximate solutions for first and second order fuctuations, respectively, using the renormalized values 
computed in section 6. We discuss the back-reaction on the geometry in section 8 and we give our conclusions in 
section 9. In the three appendices we A) compare our analytical approximation with the WKB method, B) exhibit 
the fourth order adiabatic expansion and C) compare some of our results with those obtained in a different gauge. 



II. LINEAR PERTURBATIONS IN THE UNIFORM CURVATURE GAUGE 



We consider inflation in a flat universe driven by a classical minimally coupled scalar field with a general potential 
V(<f>). The action is: 



S : 



d A xC 



d A Xy 



fJ 



R 



16ttG 



(1) 



where C is the Lagrangian density. 

Let us now study the fluctuations of the scalar field <p(t,x) around its homogeneous classical ^2 value (j>(t) and 
include metric perturbations. For the homogeneous case we have 



3iJ0 + U = 



H- 



IT 



v 



where H = a/a is the Hubble parameter and a is the scale factor. 
The scalar perturbations around a flat Robertson- Walker metric are: 

ds 2 = -(1 + 2a)* 2 - aP, i dtdx i + a 2 [8^(1 - 2ip) + 2^] dx l dx j , 



(2) 



(3) 



where the symbol ,j denotes the derivative with respect to the spatial coordinates. We choose to work in the uniform 
curvature gauge: 

ds 2 = -(1 + 2a)dt 2 - afadtdx* + a 2 5 l] dx i dx 3 . (4) 

We note that this choice fixes uniquely the gauge, just as the more frequently used longitudinal gauge (for a review of 
cosmological perturbations in this gauge see |24|). This can be seen by setting e° = —i/j/H, e = 7, where & = (e°, e ,J ) 



is an infinitesimal coordinate transformation 



longitudinal gauge we write the metric |24j in that gauge: 

ds 2 = -(1 + 2<$>)dt 2 + a 2 {l - 2y)S, lj dx l dx j 



e^). In order to see the connection with the more known 

(5) 



and note that the transformation between the two gauges can be obtained through a time reparametrization e° 
-a/3/2, e = 0: 



d 
di 



aH 



-13. 



Let us now derive the equation of motion in the uniform curvature gauge (see also [2£ 
fluctuations obey to the following equation of motion: 



(6) 

43]). The scalar field 



1 2 

ZHib -V (p + Vththip — 2aq> - 

a A 



6Ha-^-X7 2 {3 
2a 



2aV 6 



2a 



V 2 /3, 



(7) 



3 



where a dot denotes a derivative w.r.t. the time t. Starting from the Einstein equations, G% = 8ttGT^, in order to 
obtain an equation for ip only one needs the energy and momentum constraints in their linearized version, i.e. the Gq 
and G° linear equations: 



8nG 



4> 2 d (H 
~E ~dt 



-f 



(8) 



a,{ = ^G— (p., 



(9) 



Because of the absence of anisotropic stress we also have + 2H/3 — 2a /a (this relation replaces the equality $ = 
in the longitudinal gauge) [2(j. On subtituting these two latter equations in Eq. Q one obtains 



1 9 



4K 



v? = 0. 



(10) 



It is important to note that the effective potential for the fluctuations can be rewritten as: 



^ + 2§ 3tf-§ + 2 | 



3H ) 



(11) 



where we have introduced the (positive) slow-rollover parameter e = —H/H 2 . This means that the self-consistent 
inclusion of gravitational fluctuations changes the effective potential for the field fluctuations. In particular, gravita- 
tional fluctuations generally decrease the effective mass to first order in e. 

Equation (|10fl is generally seen written in conformal time. The Fourier tranform modes of v = cup satisfy: 



v'l + 



k K = , : 

z 



H 



(12) 



where ' denotes a derivative with respect to the conformal time ry, dr\ = dt/a. On comparing the last equation with 
Eq. (12) of [2j| it is immediate to see that tp satisfies the same equation as the Mukhanov variable Q. Therefore, the 
uniform curvature gauge has the advantage of singling out the true dynamical degrees of freedom (the matter ones) , 
even if it has the disadvantage of being non diagonal in the metric perturbations. 



III. BEYOND THE LINEAR ORDER 



To second order we consider a metric having the following coefficients: 

5 oo = -l-2a-2a (2) 



9o i 

m 



1(0, 

Si -, - 



(2) 



■0. 



9jX 



(2) 



,(2) 



(13) 



The above metric is the extension of the uniform curvature gauge to second order: a*- 2 -* and [3^ are the scalar 
perturbations to second order. To second order, scalar, vector and tensor perturbations do not evolve independently 
as is the case in first order. For this reason we take into account second order vector and tensor perturbations, 



represented by the divergenceless vector x) an d by the transverse and traceless tensor h\^' , respectively. In the 
above we have omitted first vector perturbations (which die away kinematically) and tensor perturbations (which 
satisfy the usual equation h + 2>Hh — V 2 h/a 2 = 0). With this approximation we are neglecting the EMT of vector 
and tensor perturbations, and their correlations with the scalar perturbations. We finally note that the choice in Eq. 
ff 1 31) (including vector and tensor metric elements to first order) fixes the gauge completely to second order. 



(2) 



The Einstein tensor expanded to second order is: 

= -3H 2 --X7 2 f3 + 6H 2 a 
a 



--V 2 /3^ + 6H 2 a^ - UH 2 a 2 + ^H 2 \Vf3\ 2 
a 4 



+ — (Va • V/3 + 2aV 2 /?) 

+^(A^-(v^) 2 ) 



G° = G° (0) +^ (1) +< (2) 

= -2Ha,i - 2HJ 2) + %Haa ti - ^-Ha tl \7 2 (3 

+ ^Va • V/3 4 - |v/3 4 • V/3 - \v 2 xf ] 
G) = Gf ] + 8Gf ] + «5Gf } 

= 5) {-(3ff 2 + 2H) + 2a(3H 2 + 2H) + 2Ha 



+J_ V 2 a--V 2 /3--Lv 2 /3 
a 1 a 2a 

+2a^(3H 2 + 2k) + 2Ha^ + -iW 2 ) - -V 2 /3^ - — V 2 /3^ 

a z a 2a 

+^Va • V/3 + y V/3 • V/3 + Q|V/3| 2 - 4a 2 ^) (W 2 + 2tf) - 
-8ffad + f A + V 2 /? - 4"V 2 a + -V 2 /3 - -^|Va| 2 



2a 



+ ^ ( ^ A ,„_ (vW) } + {^, + ^,_J_„, 

2a ' J a ,J a" 2 ,J a^ a ar 



1 



2. (2)i 



To second order the scalar field is expanded as: 

0(t,x) = <^(t)+^(t,x)+^ 2 )(t,x) 

The EMT of inflaton fluctuations 
to second order is: 



1 o 



= T 



0(0) +5T 0(l) +(5T 0(2) 



1 ■ 



-J 1 + v{<t>) 



+ <p a - (flip - V^ip 



+0V 2 ) - fom - 1>( 2 ) - 1^ 2 - J-^i 2 - i W 



-20 2 a 2 + 2(j>aip + 



1 V 



|V/3| 2 



5 



rpl 



(2) 

T m +S T^ + ST^ 



1 V 



jxp-<f> a - V^ip 
1 • 



&p< 2 > - 0V 2 > - l/ ^ 2 > 



+i f <yi 2 - 4^0 - 7<^ 2 |V/3| 2 + 40 2 a 2 + -<j>Vtp • V/3 - -^| V^| 2 - V,^ <p' 
2 \ 4 a 



(20) 



(21) 



IV. EQUATIONS TO SECOND ORDER 

In this section we exhibit the second order equations corresponding to Eqs. I|8I9I1U|) . using, when convenient, the 
homogeneous and first order equations. We first give the momentum constraint to second order, which is obtained 
from the expressions given in the Eqs. (|15l) and (|20() : 



(2) 

where \i are vector metric perturbations and 

4ttG 



s, 



s,= 



H 



(ip - 20 a) Vti + 4aa,i - (a *V 2 /3 - fly) - lp* ft, 



The term Si can be written as + Vi, further 



(22) 



(23) 



= AnG—ip 
H 



(2) 



where 



-47rGe^ 2 + 2a 2 



;|V/3| 5 



1 

"V2 



AttG. 



H 



AaH 



(24) 



(25) 



contains the quadratic contribution of first order perturbations. We note that s includes non-local spatial contributions 
which nevertheless to leading order in e and for long wavelenght have an ordinary behaviour on large scales. One can in 
fact approximate in such a limit, for the isotropic case, the first term in the second line of eq. (|25|l as (2irG) / H (pip. On 
combining Eqs. (|14I19|I . the energy constraint to second order, upon using Eq. I|24|l and the lower order constraints, 
is given by: 



a H dt 



H 



,(2) 



Q + UirGVs , 



(26) 



where Q is defined as: 

Q = 12ff 2 a 2 -^ 2 |V/3| 2 -^ 
-8?rG 



\p 2 + ^2 I v ^| 2 + \ V ^ 2 + 2 ^« 2 " - V/3| 2 



(27) 



We observe that Eq. (|26|) for (3^ is reminiscent of a universe filled with two components, with the term Q — IQttGVs 
playing the role of a non-adiabatic pressure term |24j |. upon comparing with the first order equation (JSJ. From this 
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analogy one may guess that the terms on the right hand side of Eq. I|26|l should approximatively cancel. Explicit 
calculations confirm this property. 

The equation of motion for the scalar field to second order, after using all the previous constraints, is given by: 



(2) 



J_ vV 2) 

a 1 



4K l-I 



<pV=D, 



(28) 



with an homogeneous part which is the same as found, see Eq. (|1(JII . for the first order fluctuations and an inhomo- 
geneous term D which is purely quadratic in terms of the first order fluctuations. In particular one obtains 



and 



R 



D = R + <j)s- 2V 4 ,s + 7-77 (Q - 16nGVs) 
2H 



1 2 cb 1 

V^ U Lp 2 -2<j)aa + aip+ -»aV 2 ip - 2V^aip + ^-a ti (3' 1 - -V^ (3\ 2 + 



(29) 



2 
1 



a 



2a 4 



2a 

2a a 



(30) 



We have therefore obtained the perturbative equation for the scalar field fluctuations to second order, which is a 
novel result. 



V. APPROXIMATE SOLUTION FOR LINEAR PERTURBATION 



The equation for the Mukhanov variable (|f 2f> does not have an exact solution, except for the known case of an 
exponential potential [26|. Approximate schemes to obtain the long wavelength solution, such as the slow- rollover 
technique |27j . exist. Any approximation must agree with the solution to Eq. (jl2[) with k = 0: 



/dt 
2 
az z 



(31) 



where C and D are constants (when the first term in Eq. I|31|l dominates, the curvature perturbation £ = Hip/4> is 
constant to leading order). For small wavelengths (k ~ 0) C and D include the k dependence of the modes. However, 
in order to perform an adiabatic regularization, we do not only need the solution for small k, but for the whole 
spectrum. The slow rollover approximation applied to Eq. I)12|l leads to a Bessel function whose infrared limit agrees 
with Eq. H31|l only on freezing a spurious time dependence of the Bessel solution to the value on horizon crossing. 
This problem is the same as was encountered in 17J on trying to extend the approximate solution for the inflaton 
fluctuations in a rigid space-time to fields with different mass. Therefore, the slow rollover paradigm is not useful in 
order to obtain an approximate solution over the whole k spectrum. 

In order to consider the slow-rollover of massive chaotic inflation we take H ~ — m 2 /3, H ~ and neglect terms 
0(H 2 1 H 4 ) ^tJ- During slow-rollover the equation for ip is then: 



ip k + 3Hip k + 



m, 



6H 



•Pk — . 



(32) 



On comparing with the fluctuations in a rigid spacetime considered in [ljj , we see that the gauge- invariant fluctuations 



in Eq. (|1(J|) have a negative mass since H 



2 /3: this is due to the proper inclusion of gravitational fluctuations. 



This result is true for many inflationary models and it is related to the red spectrum of the curvature perturbations. 

In order to have an approximate solution for the whole k spectrum we proceed in analogy with [17] and we choose 
for large k: 



1 



2 3/2 



7rA 

AH 



1/2 



(33) 



where £ = k/(aH) and 



A = 1 — 
3 

" = 2" 



H 
1 m 2 

zip 



= 1 



H 



2 m 2 

zip 



2e. 



(34) 
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We note that v for the gauge-invariant fluctuations differs from the corresponding quantity for inflaton fluctuations 
in rigid space-times [l7j . 

The reason the procedure followed in jljj does not succeed in producing an approximate solution, valid over the 
whole range of k, is the dependence of v on H. In fact, among the terms which are apparently 0(H 2 / H 4 ), the 
following term: 



leads to a term which is of order 0(H/H 2 



\ k \(H 2 -H 2 



(35) 



(36) 



making the approximation not valid for small k and at large times. For this reason we consider Eq. i|33|) only as an 
approximation for modes which are still inside the Hubble radius. We also note that the solution for the ultraviolet 
regime m Eq. (J23l is the same as that of the slow-rollover approximation, although here H ~ const, and not e ~ const. 

On considering the infrared limit for the gauge-invariant fluctuations in Eq. (|31J) we know that tpk ~ 0(1/ H), since 
(j) ~ const. For k < aH we must to replace the solution in Eq. l|3*3*j) in order to reproduce the correct behaviour in the 
infrared by: 



,3/2 



7rA 



1/2 



H(t k ) 
H(t) 



#3/2 ( A ' 



where tk is a time related to the istant for which the mode k crosses the Hubble radius, 
estimating the time tk to appendix A and give the Hubble parameter for this time value 



(37) 



We leave the details of 



H(t k ) = H a 



2 m log 



(38) 



We note that we have a nearly scale invariant spectrum due to the dependence on tk- 

We also note that the amplitude of the curvature perturbations, = Htpk/(p, associated with the solution in Eq. 
()37J) . need not obey the observational constraints (fc 3 ^ 2 |Cfe| ~ 10~ 5 ). The reason is that the spectrum is red tilted 
and we are considering modes which could have exited the Hubble radius much earlier than the ones relevant for 
observations. Only if the duration of inflation is minimal in order to solve the horizon problem do the curvature 
perturbations associated with Eq. (|37(l satisfy the observational constraints. 

Cosmological fluctuations are canonically quantized as usual and we shall consider the vacuum defined by Eqs. l-i-il) 
and l|37j> : 



<£(t,x) 



1 



(27T) 



3 

k 



[pk[t)e 



ik-x 



-ik-x 



(39) 



where the bk are time-independent Heisenberg operators. 



VI. THE ADIABATIC SUBTRACTION 

We now compute the integrals obtained on taking the expectation values of the relevant operators. The strategy 
is slightly different from that of our previous calculation jljj since the mode function is obtained by matching two 
approximate solutions. Therefore we split the integral: 



r+ca 


v aH 




/ dk -> 

h 


[ A dk + 

J I 


l 



dk 



(40) 



where t — CHq is an infrared cut-off related to the beginning of inflation 44] 28, 29] and k = is the turning-point 
of the Bessel function in Ea. (|33|l . We substitute the infrared solution (|37() in the first integral and the ultraviolet 
solution H33JI in the second integral. For the ultraviolet part we employ, as in the previous article 0]) dimensional 
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regularization using a d dimensional space measure. Subsequently an adiabatic subtraction is performed in order to 
obtain the renormalized quantities. 

From our previous calculation we know that the leading contributions come from terms which contain (<p 2 ), while 
terms such as (V^tpV^tp) are more ultraviolet and therefore non-leading. By the symbol (...) we denote the average 
over the quantum state defined by Eqs. and l|37|) . It is important to note that the standard fourth order adiabatic 
subtraction is sufficient to regularize the EMT of cosmological fluctuations as well. In fact, the terms which apparently 
would not be regularized by a fourth order expansion (see the term ^ — (V 2 /3) in the Gq equation i|14|) . for 
instance) vanish on averaging over an homogeneous state. As an example of the calculations we exhibit the details of 
the calculation of (^)ren- 

The ultraviolet and infrared integrals are respectively: 

, imiv R 2 f „ ™ 2 ^H 2 ( 2m 2 \ / 1 \ d ~ 3 { aH\ d ~ 3 (\ d\\ 

^ = m? H \ A + b Jp + c ^ - ( 2 + 3 W) ("x) r U ~ 2) ) + 

+0(d-3), (41) 




where eo = ~H/Hq an( i A, B and C are constants with a complicated dependence on hypergeometric functions and 
we have taken for the calculation of Ea. (|42|) £° ~ as is done for the case of the t/. derivation in appendix A. At 
the end of inflation log a — > l/(2eo) which is the inverse of the slow-rollover perturbation parameter. Therefore in the 
above expression, in order to be as accurate as possible towards the end of inflation we have included all the leading 
log terms which are multiplied by the small parameter eo- The adiabatic fourth order is (see appendix B and |Tt| ): 

2 H 2 J 127 43 to 2 29 H 2 ( 2m 2w 1 x d ~ 3 

W )(4) - I -^--^IP + 1W - [ 2+ 3H~ 2 J V2tt 1 /2 / 

(am) d - 3 T Q - + O [^j } + 0(d - 3) , (43) 
and the resulting renormalized quantity is 

= (M' H2 ^ { i 1 - 2f " log ^) 3 - (1 - 2e ° log ° )3 } + 

ft rri f a r, m2 ^ 2 J 1 rn 2 \ , ( vH \ 127 43 m 2 29 H 2 

H 7tH i A + B^r + C^T - 4 H r k 1 „ T + 

16tt 2 \ H 2 m 2 V ZH 2 ) \\m) 45 90 H 2 15 m 2 

The leading behaviour of the renormalized correlator at the beginning of inflation, which actually originates in the 
infrared region, grows with log a: 

(^W-^W^Vlna. (45) 



Its time derivative is 



4tt 2 V H , 

d , 9 , ft Hn , 9, ft Hn 2 m 2 , 9 , , , 

jyUv. ~ ^ + 2eH(^ = + --(^ 2 ) REN . (46) 

We note that this result is very different from the de Sitter result (H constant in time and e = 0) since the second 
term on the right hand side may dominate for large times. 

The above result in Eq. I|45|l is reminiscent [45j of the usual Hartree-Fock term formulae for stochastic inflation 
(see Eq. (5) of [HI): 

5<^— j? *-^— (47) 
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However, Eqs. i|46|) and 1(47(1 differ not only in the driving terms, but in particular in the mass term: in our result 
there is a negative mass term, whereas in the stochastic approach in rigid de Sitter space time there is a positive mass 
term. Our result 1)46(1 is completely consistent with Eq. 1)32(1 where the gravitational terms change the sign of the mass 
term in the evolution equation for the fluctuations. The difference is therefore due to the inclusion of gravitational 
fluctuations which are not properly taken into account in stochastic inflation. 



At the end of inflation, for the characteristic time scale of slow-rollover St 
in l(44|) saturates to 



3Ho/rn 2 , the correlator in the expression 



2\MAX 
/REN 



47T 2 



Ho 
H 



1 

6e 



h H$ 1 
8^ H 2 m 2 



(48) 



It is instructive to also compute (</?</?)ren and (^ 2 )ren by the same method, using the approximate relation for the 
infrared modes tp k ~ —(H/H)tfk- The final result (omitting all the steps in Eqs. 1(411421431 144(1 ) is: 



H((pip) KEN = H(ipipy H + 



16tt 2 



:H 4 



' ^m 2 8 m 2 , ( vR 

D + E— + In 

H 2 3 H 2 I Am 



4 m 2 / 1 
9H 2 + W 



(49) 



REN 



2\IR 



16tH 



:H 4 



„ „m m , / vH 



61 27 m 2 ( I 

60 ~ Toll 2 + 



(50) 



where D } E, F and G may be written in terms of generalized hypergeometric functions and for the infrared contribution 
one has |46| 



IR 



He(<p 



2\IR 



(51) 



2\IR 



2,2 l,J2\IR 



(52) 



VII. APPROXIMATE SOLUTION FOR SECOND ORDER SCALAR PERTURBATIONS 

An approximate solution for the leading quantities to second order can be obtained for large scales where the 
infrared modes dominate. Let us consider ip^ 2 ' which satisfies Eq. ((28(1 . On neglecting terms 0(H / H 3 , H 2 / H 4 ) and 
the second order derivative we obtain: 



ZHtpW _ m 2 ^(2) 



V, 



2 2 

in tp z 



(53) 



which can be rewritten, on also using the homogeneous equation of motion for the inflaton, as: 

d 



3 ft K 2) ) * * 



2 2 

2HM* 



(54) 



Let us note that we use the notation M 2 X = l/(8irG) for the (reduced) Planck mass definition. In the same approxi- 
mation the energy constraint is: 



2a^(3H 2 + H) ~ -6^ " 2 
'pi 



-fow - - l -v HV 2 



(55) 



The expression we want is (p^), which is averaged over the initial vacuum state, and can be obtained by using 
the leading contribution in Eq. 14511 to the quadratic correlator of the first order fluctuations (^ 2 )ren- We obtain 
directly from the last equation: 



(2)\ 



h -Him 2 



1 1 H_ 

AH 3 + 2HH 2 n M 



(56) 



where M is a constant with the dimensions of a mass. It is also convenient to exhibit a slightly different form, which 
is of course equivalent in the large a limit at the end of inflation since the infrared contributions are then a maximum. 
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Since <j> is almost constant one may also obtain, from Eq. (J52J, 

/ „(2)\MAX ^ 9 ;,„2\MAX /c 7 \ 

It is also useful to write the contribution to (a^ 2 ^) which we shall later use to study the backreaction effects on some 
scalar observables. One obtains for its leading behaviour 

/„(2)\MAX^ Y / (2)\MAX 3 1 ^ /.JlMAX ^ e / 2\MAX /ro\ 

{a } ~2HM^ X - i ^^(^)REN-^(^>REN (58) 



VIII. APPROACHES TO THE BACK-REACTION 



One may follow different approaches in order to study the backreaction effects due to cosmological fluctuations. In 
the following we shall consider two methods in order to tackle this issue. 

One consists of considering only the first order perturbations, then imposing to first order the energy and momentum 
constraints and finally defining an effective EMT by including all the quadratic terms present in the Einstein equations. 
Subsequently one averages everything over the quantum vacuum and first order quantities disappear. One finds that 
the effective EMT which appears in the averaged Einstein equations is modified by the back-reaction. The results 
obtained in sections [HI and IVII are sufficient for this purpose. However one still has to study how any observable 
averaged over the quantum vacuum is then affected. 

A second approach is related to the standard perturbation analysis of the Einstein equations up to second order. In 
this case we impose the energy and momentum constraints and study the inflaton equation of motion perturbatively 
up to second order. One does not define any modified EMT but directly studies in this framework any observable 
averaged over the quantum vacuum. In order to pursue this approach results given in sections IIIII IIVI and IVIII are 
also necessary. 



A. The Energy-Momentum Tensor of Cosmological Fluctuations 

By the EMT of cosmological fluctuations we mean the second order part both of the scalar field EMT and of the 
Einstein tensor which is quadratic in first order fluctuations: 

__ rjif_t quadratic j\/f^ ^f/i quadratic (59) 

This method of considering the EMT of gravitational fluctuations is treated in textbooks |32j and has also been 
previously used in In this scheme one considers the modified Einstein equations M^G^^ = + (rjj) 

which therefore include back-reaction effects. 

Let us consider the EMT of cosmological fluctuations averaged over the state annihilcd by the operator b defined 
in Eq. 

For a generic potential the leading terms in the energy density are 



It is important to note that the second term on the right hand side is in general positive and larger than the first 
during inflation. This second term is the contribution of metric perturbations. On using the slow-rollover parameters 
e and r\ (= M 2 ,-^-) we can rewrite the leading terms in Eq. (|60|l during slow-rollover as 

(r }^-e~-^V)(l-4£) . (61) 

Analogously the average pressure is: 

(rj) P S) - 6] (-^ 2 ) + 12HHll x {a 2 )) ~-eS) (62) 
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We now restrict the analysis to a quadratic potential with 77 ~ e. On using Eq. i|45|l the leading behaviour for the 
initial time of the energy and pressure density is: 

£REN ~ -PREN ~ -^ 2 (V 2 >REN ~ --^m 2 H 2 (jjf^J m a = ~ ^~2 eH ln G ' ( 63 ) 

further, on using Eq. l|48|) . we have the following maximum value: 

max max 3 ^ H o (f - A \ 

e REN ~ "PREN ~ Jp 

we then obtain that for 

H ~ (16^ 2 /^) 1/6 ^ 2/3 ^ P Y 3 ■ (65) 

the EMT tensor of cosmological fluctuations cannot be neglected with respect to the energy content of the homogeneous 
classical inflaton. 

As a final remark we note that the amplitude and time dependence of the EMT depends on the approximation 
used for linear perturbations. The slightly different results in the amplitude of the EMT with respect to 0, are 
explained in Appendix C. 



B. Back-Reaction on the Geometry 

In the perturbative approach to the Einstein equations any back-reaction effect is analyzed by evaluating pertur- 
batively quantities which characterize the geometry. 
The expansion scalar is defined as: 

6 = V„u" , (66) 

where is a four- vector field defining the comoving frame. The four- vector is normalized, u^u^ = — 1. To second 
order the four vector is: 

1 ^ 3 2 (2) i P- l + P {2) •* (f . 7 , 

u = 1 — a H — a — cr ' , u — (67) 

2 a 

and the expansion scalar reads: 

= 3H- 3Ha + -\7 2 fi - ZHa {2) + -V 2 (i {2) - 
a a 

+ ^Ha 2 + -Va- V/3 + • V/3 . (68) 
2 a 4 

The leading terms, once we average, are: 

<e) = e(°> + <e( 2 >) 

= 3 H (l+ 3 -(a 2 )~(aW) 

(69) 

where the results given in Eqs. 1451) and l|58|l have been used, together with the relations a = <\>j (2HM 2 l )(p and 

(f> 2 = —2HM^ l , in order to calculate the only two non negligible second order contributions. We note that by 
leading terms we mean the first corrections in e which are also leading in log a, thus essentially including the infrared 
contribution present in the renormalized quantities, as already discussed. We observe that the result in Eq. (|69|l is 
completely consistent with the considerations on the EMT made in the previous subsection. When the magnitude of 
the EMT of cosmological fluctuations is of the same order as that of the background (0) changes significantly. 
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The result in Eq. is in contrast with where a vanishing result was obtained for the back-reaction on the 
expansion rate (0) and where only first order classical perturbations are considered. We note that also in our gauge 
choice 8 is related to the Einstein tensor for long wavelenghts and leading order in e: 



6 ~ 3F ( 1 - a - a {2) + ^a 2 ) ~ yj -3G° , (70) 



and therefore the result in Eq. 1|69|) could also be obtained by expanding y/— 24ttGTq to second order and averaging 
over the quantum vacuum. The difference with respect to [l2j may be due to the absence of second order fluctuations 
in ^3] and/or to a peculiarity of inflation driven by a quadratic potential. 

We now consider another geometric scalar observable which is associated with the rate of change of the expansion 
scalar, i.e.: 

n = t^VvV^u" . (71) 

Its expansion up to second order is given by 

n = n^ + ci^ + n^ 

H 1 

= 3H -6Ha-3Ha V 2 /3+-V 2 /3 + 

a a 

-6Ha^ - 3Ha® _ ^ V 2 /3 (2) + ±V 2 /3^ 
a a 

H H 
+12Ha 2 + \2Eaa - 4— Va ■ V/3 + — aV 2 /3 - 
a a 

1 o • 1 d -j A \ Id 



i^ + i|(^.w) + i|(v«.v ? ) 



+ \^fi ■ V (V 2 /3) (72) 

For the homogeneous case G = 3H and therefore £1 — 3H. During slow-rollover for a massive inflaton H (~ — m 2 /3) 
is constant in time, and therefore is a gauge- invariant quantity up to second order [|| l33| . 
The leading terms, once we average, are: 

(fi) = n<°) + (n( s ') 

= 3H- 6H(a (2) ) - 3H(a {2) ) + l2H(a 2 ) + 12H{aa) 

~ 3fr(l + F(e,logo)) (73) 

where F(e, log a) ~ eF\(H) + e 2 log a F2(H) + Therefore one observes that in this case the leading corrections 
to order e which come from the infrared region cancel, while these were present in (O). Let us again note that at 
the end of inflation driven by a massive inflaton log a ~ l/(2eo) so that any logo correction effectively reduces the 
perturbative order in powers of e of the correction. 

The result on (fi) should be interpreted with care since we are not able to go to higher order self-consistently 
within our approximation. We also note that this vanishing result may also be a peculiarity of inflation driven by a 
quadratic potential. In any case the rate of change of the Hubble parameter can also be studied more carefully at the 
background level, without considering any back-reaction. Including the first non trivial corrections one finds 



H^-^ l-l%). (74) 

To conclude this section we stress that in both the results (|63|) and (|69|l the corrections to the background values, 
which grow with time, are 0(e). In the de Sitter limit (e = 0) these corrections vanish, consistently with the absence 
of physical scalar perturbations to first order for a universe driven by a cosmological constant. 



IX. DISCUSSION AND CONCLUSIONS 



The renormalized EMT of cosmological fluctuations during inflation with a quadratic potential is studied. With 
respect to our previous work |17| we have self-consistently taken into account the gravitational fluctuations accompa- 
nying the scalar field fluctuations. 
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We find that the renormalized EMT of cosmological fluctuations during slow-rollover carries negative energy density 
and has a de Sitter like equation of state to leading order. The negative sign for the renormalized energy density is 
due to the inclusion of gravitational fluctuations (see also 0, El for the same claim regarding long- wavelength modes) 
and did not appear in the previous calculation in rigid space-time This effect is generally true for single field 

inflationary models with e > rj/4, as is seen from Eq. I|61|) . 

The back-reaction problem was also treated up to second-order in perturbation theory. For this purpose we extended 
the UCG gauge beyond the linear order by taking into account vector and tensor perturbations to second order. The 
inclusion of vector and tensor perturbations to second order is important in order to obtain the correct equations 
for the second order scalar perturbations. We derived the equation for the second order scalar field fluctuations 
and approximatively solved the system of second order equations for single field inflation during slow-rollover. We 
computed the expansion rate and its coordinate invariant derivative f2 to second order, and we found that the 
expansion rate is decreased by the EMT of cosmological fluctuations, while fl is not affected to leading order. A non 
vanishing correction for f2 will probably appear in the next order. 

We do not find evidence that the back- reaction of cosmological fluctuations is a gauge artifact as is claimed in ^jj • 
The second order corrections to the local expansion rate are non vanishing and are computed to lowest order in Eq. 
(16911 . We note that the back-reaction of cosmological scalar fluctuations vanishes in the de Sitter limit (e = 0), as it 
should. Since is not a gauge-invariant quantity in the model under consideration [ill8.'i|. we also computed in the 
gauge wherein the inflaton is a clock: such a check was done by a completely new calculation (illustrated in appendix 
C) in order to also study the transformation properties of the metric elements. We confirm that back-reaction is an 
0(e) effect, but its time behavior is different (there is no log a enhancement from the infrared region) because of the 
differing times in the two frames. In this new gauge calculations are also more involved since a higher accuracy in the 
solution for the Mukhanov variable would be required. 

It is interesting to note that the inclusion of gravitational fluctuations sistematically contributes negatively. To 
first order gravitational fluctuations contribute negatively to the effective mass for the fluctuations, and to second 
order carry negative energy density. The self-consistent inclusion of gravity also changes the stochastic picture in an 
inflationary background. 

The approach we have used is mean field theory on a curved space-time. The procedure of averaging over a quantum 
state with a cut-off related to the patch which undergoes inflation can be seen as a spatial average over the particle 
horizon. This avera ge i s completely different from the average over the Hubble region during inflation used in the 
stochastic approach [3(| and in which gravitational fluctuations were neglected. We think that the difference between 
the results of stochastic inflation and those found here are not just due to a different coarse graining, but are also 
due to our inclusion of gravitational fluctuations. We think that the relation between these two different approaches 
should be further investigated. 



We wish to thank A. A. Starobinsky for discussions and suggestions on the draft. We also wish to thank R. Abramo, 
R. Brandenberger, A. Linde and S. Matarrese for comments. F. F. thanks R. Abramo and R. Brandenberger for many 
conversations on the back-reaction problem over the previous years and N. Bartolo and B. Losic for discussions. 



In this appendix we show how the procedure used in this paper to match solutions which are good approximations 
in different ranges of k space includes the recent application of the WKB method in the context of cosmological 
perturbations [34l l35| . In this paper we have chosen a Bessel form in Eq. (|33l) for modes which are still inside the 
Hubble radius. The scaled equation [ip^ = a?/ 2 ipk) associated with eq. <|32[1 is given by 
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X. APPENDIX A: COMPARISON WITH WKB METHODS 




(75) 



where v 1 = v j\. The turning-point of this equation is obtained for 



a(t) 



k 



(76) 



H(ty(t) ■ 
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So, for our calculation, we can neglect the weak time dependence on the right side of Ea. l76|) with respect the 
exponential dependence on the left side and obtain the following result for the turning-point 



tk 



Hn 



1-4/1 + 2 



H 



1 X Q k 
— log 

Hq HoVq 



Ho \ V Hq - Hova 
The solution near the turning-point of Eq. i|32[l . on using the WKB method, is given by 



1 



Si At ( al /3 (t - t k )) + B 2 Bi (al /3 {t - t k ) 



2 3/2 

Around the turning-point we can consider the following approximate expression 36] in Eq. 



H<P (x) = ^ 



( w— —z — -arctan w I 



oil 

V 



(77) 



(78) 



(79) 



S« (a 



S arctan toy j ^ 



tr- 



im 



v\ 3 



w 



(80) 



where w — 'x 2 — v 2 . We have checked numerically that the above approximations are remarkably good. The Airy 
functions are precisely the functions used to match the WKB solution, ea. (j78|l . inside and outside the Hubble radius 
|35|. In particular around the turning-point the match is specified by 



l/ ( t k)\" 2_ 1/3,. , , 

1 w « a k ' [t-tk) 



(81) 



Therefore, by choosing a Bessel form inside the Hubble radius one already has the correct form in the vicinity of 
the turning point. This cannot be achieved with the WKB approximation, since the WKB solution is singular in 
correspondence of the turning point. 

XI. APPENDIX B: THE FOURTH ORDER ADIABATIC EXPANSION 

In order to remove the divergences which appear in the integrated quantities as poles in the T functions, we shall 
employ the method of adiabatic subtraction [33 • Such a method consists of replacing our functions with an expansion 
in powers of derivatives of the logarithm of the scale factor. This expansion coincides with the adiabatic expansion 
introduced by Lewis in [38] for a time dependent oscillator. 

Usually it is more convenient to formulate the adiabatic expansion by using the modulus of mode functions Xk — 
\<Pk/V2\ and the conformal time r\ |37| (dr/ = dt/a). We follow this procedure and write an expansion in derivatives 
with respect to the conformal time (denoted by ') for Xk- We then go back to cosmic time and insert the expansion 
in the expectation values we wish to compute. Adiabatic expansions in cosmic time and conformal time lead to 
equivalent results, because of the explicit covariance under time reparametrization j3fj. 

The variable x k satisfies the following Pinney equation: 



Xk + 3Hx k + 



rfc 2 



m 2 + 6H 



1 



Xk 



a 6 xl 



Following p"?! we rewrite Eq. I|82() in conformal time in the following way: 

1 



(ax k )" + £l k (ax k ) 



(axkY 



(82) 



(83) 



where 



2 7 2 1 2 2 1 2n 

— k + ma an 

6 



(84) 
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and R is: 



R = R — 36H . 



(85) 



with R = 6-\ the Ricci curvature. From Eqs. (|83I84I85|) one obtains the expansion for Xk up to the fourth adiabatic 
order: 



1 1 



' ^ V 1_ 4 e2+ 32 e2 "4 £4 



(86) 



where f2fc is defined in Eq. I|84|) and £2 , £4 are given by: 



(■2 



1 fi^ 3 fi' fe 2 
'2Q| + 4 



1 » 

64 " 4fif 2 4r»f 2 



(87) 



The solution in Eq. I|86|) must be expanded again since R is of adiabatic order 2. Therefore £C fe is 



1 If 1 £ 1 5 2 ^ 2 1 

^(4) " c i/ 2 E i/ 2 I 1 + I C "6 Ef + 32 C 36 Sf 



1 1 

"16 E 4 

5 1 

64 E£ 



I? 

g" 



2c 



i? 



;_ 6" 



'2 4 
C TO 



2c TO 2c TO c — 

6 6 



9 1 R „ 2 
+ 64El C 6 C m 



65 1 i? , 2 4 u 2 
256 EfV m + 32 62 



5 2 1 1 
32 ' 2 * 4 £4 *J 



where c 



and 



St 



/; 2 1 2 2\l/2 

(k + am) 1 



IE, 3 E 



'2 



£2* — — 



£4* 



2^1 
4E3/ 2 



1 1 n 

is! 62 



(89) 



XII. APPENDIX C: COMPARISON WITH CALCULATIONS IN DIFFERENT GAUGES 

a. In 0, 0] a different procedure was followed. The calculation was performed in the Newtonian gauge (see Eq. 
(O)and the regularization of the metric perturbation ($ 2 ) was considered. We can compare our results by using the 
relation between field fluctuations ip in the uniform curvature gauge and field fluctuations tp N in the Newtonian gauge: 



Assuming that <£> in negligible on large scales, we obtain |24j: 

'e- 1 



which, on using the approximation used in [7J, |9j, |24j , leads to 

e- T 



fc~ 3 / 2 



2y/6M p i 



1 - 



\og[k/(aH)) 



1 + log[a end /a] 



(90) 



(91) 



(92) 
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where a en d is the scale factor at the end of inflation. This approximation for the gauge-invariant field fluctuations is 
different from our Eq. (|37[) . 



b. After this paper appeared in the electronic arXiv, it was suggested that corrections to < > should disappear 
in the gauge wherein the inflaton is a clock, homogeneous in space. Let us note that the inflaton can work as a clock 
only when it is not oscillating, as is the case in the slow-rollover regime. 

In order to check this claim we performed a new calculation, since the gauge where the inflaton is homogeneous is 
different to the one used in this paper. The gauge in which the inflaton is unperturbed is described by a metric with 
three scalar degrees of freedom and, again, this is a choice which fixes the gauge completely because of the presence 
of the scalar field. As before we do not pay attention to the vector and tensor degrees of freedom and therefore we 
consider the following second order line element depending only on three scalars: 



ds z = -(1 + 2& + 2& {2y )dP - a{(3. 



%(l-2^-2Vi( 2 )) 



dx % dx 3 



(93) 



where we have taken a(t) = a(t), since in the perturbative approach a(t) and d(t) satisfy the same equation in t and 
t, respectively (however the scale factor is not observable for the case under consideration). We expect that ip will 
satisfy a dynamical equation of motion in this gauge. 

Instead of solving the system in this gauge, we find solutions for the metric perturbations in (|93|l by using its gauge 
relation with the metric in (|13[l . An infinitesimal coordinate trasformation up to second order |6j: 



= (i) 



£ (l),K e (l) + £ (2) 



(94) 



(where em and e(2) are the coordinate changes to first an second order, respectively) induces the most generic change 
in a geometric object T = + + T^: 



j-(i) _> fW _ j<(i) _ £ £(i) T (0) 



(95) 



T {2 



rf(2) _ rp(2) _ n rp(l) ,\(r2 m r rp 

1 1 '-Hi) 1 + 2 V e (i) e M c 



(96) 



The time reparametrization which relates the two gauges can be found by imposing that the field perturbation is zero 
to first and second order in the metric (|93|l . According to Eqs. (1951961) the field transforms as: 



ip -> Lp = tp- 



(97) 



(2) 



£ (D l e (D^ 



(98) 



which leads to 



: (2) 



5£ = ! 

6 H 



(2) 



CC 1 H 
H 2 2 HH 2 



c 2 , 



(99) 
(100) 



where we have used the gauge invariant curvature perturbation £ = H(p/(f>, which is constant on large scales and the 
expression in eq. 1)57(1. From the above we understand that we need a large time reparametrization in order to keep 
the scalar field homogeneous in space. 

The above trasformation, up to second order, is not the one required since it would lead to the presence of a scalar 
second order contribution in the traceless part of the metric which depends on a quadratic form in di(3 and diC^y 
In order to kill this contribution a scalar second order transformation in the spatial part of the coordinates is also 
required, e| 2 % = d l e s , 2 \, while still having e!^ = 0. We find 



+ - 



3 d l d J 



2(V 2 ) 2 



-2^(D D 



y £ (l) 



(101) 
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where D^j = didj — 1/3V 2 <%. 

On using the gauge transformation i|95|) we have the following relation for the metric fluctuations to first order: 

2 



a = a-e° {1) , /3 = - ~e° (1) , ^ = He\ iy 

On solving for the metric fluctuations in the UCG gauge, using Eqs. (18191 we have 

V 2 /3 



whereas in the unperturbed scalar field gauge we have: 

C V 2 /3 



a 



H 



v 2 c 



(102) 



(103) 



(104) 



From the last relations we see that the cost of keeping the scalar field homogeneous in space is to squeeze a large 
fluctuation in the metric, i.e. ip. In the UCG gauge metric fluctuations are suppressed with respect to i.e. 0(eC , e C)j 
whereas in the gauge l|93[) they are not. To second order we have: 



« (2) = a( 2 )-e° 1) d-2ae° 1) -% + k ) et 1) +e^ 



;-0 



= a< 2 > 



(H) + IA c2 + il + ICC 

2 2ff 3 ^ H 2 2 H 2 



where, from eq. (|58|l . one has a^ 2 ^ ~ 2e 2 £ 2 . Moreover we find for [3^ a transformation such that: 



a £ (2) + ae (2)+ 2adt ( e (l)) + 



& 



(105) 



(106) 



Let us remember that we chose not to exhibit in detail the vector and tensor degrees of freedom, but it is straightforward 
to do so. It is suffices to say that to kill the possible vector degrees of freedom arising in goi, a non vanishing vector 
component in the coordinate tranformation to second order, e^j, has to be introduced, but this does not affect the 

scalars to second order, a*- 2 ), (3^ and ^ 2 \ 
Finally, for this last scalar we find 



2 s 
£ (2) 



(107) 



We are now ready to compute the expansion rate O, which is a gauge dependent observable, in the unperturbed 
(uniform) field gauge (UFG), while neglecting vector and tensor contributions and only keeping terms up to second 
order: 



= 3H- 3H&+ -\7 2 f3- 3-0' - 3Ha {2) + ^V 2 /3 (2) - 3^ (2)/ 
a a 

+ -Ha 2 + 3ot// - 6^' + -Va • V/3 - -V/3 • + -V(3 • V/3' . 
2 a a 4 



(108) 



Let us observe that the prime denotes the derivative with respect to t, and the relation between the two time 
derivatives may be written as 



d_ 

dt 



~2 di£ ( 2 )fe~- 



_I e 2 1 17 



1 d 3 

1h 2 H 2 + AHH 2 dfi 



1 H 



H 2HH 2 



cc 



1 

W 2 ' 



d_ 

dt 



(109) 



Clearly the two time derivatives coincide on acting on second order terms if one does not go beyond such an approx- 
imation, however one must exercise care for terms of lower order on going from one description to the other. 



18 



Further the difference beetwen the spatial derivatives is given by 



d 

dx l 



d 

dx l 



d_ 



(110) 



which, however, is not important since we shall neglect spatial derivatives. 

One may wish to write the corrections to the unperturbed O = 3H in the x frame, while expressing the corrections 
in terms of the already known quantities (in the x frame). For this case one must be carefull in analyzing to second 
order the fourth term of the above expression, —3dip/dt. 

To leading order, on averaging and neglecting terms with spatial derivatives, we have: 



3H 



3H 



3H 



1 



H 



(CO 



H 2 



<c 2 



-i-(5V>) 

H 

1 H 

~ 4lP 



H 



ti 



<C 2 



-(w)- 7 -( ( p 2 )+0( e 2 ^ 2 )) 



H 



H 



3H 



2M: 



(111) 



where on going from the first to the second line the terms O(CC) cancel. The infrared contributions to the terms in 
the last two lines vanish, leading to corrections to the expansion rate in the UFG gauge of order 0(e), which are of 
the same order as the neglected gradient terms. The corrections to the expansion rate in the two frames are therefore 
of the same order in the slow-roll parameters. 

If we choose to rewrite the averaged expression for 8 in terms of t, starting from Eq. (|108fl 43 an d using the 
relation in Eq. (|lt)9(l (omitting negligible spatial derivatives) to also compute H(t) = H(t)dt/dt, one finds 



(6) = 3H 1 - e 




AMI 



(112) 



which is completely consistent with our result in Eq. (|69|) and with the fact that is a scalar (but not gauge invariant), 
which, therefore, transforms according to: 



Q(x^ 



(113) 



We have then checked that the back-reaction of scalar cosmological fluctuations is 0(e) in the expansion rate: such 
a correction indeed vanishes in the de Sitter limit, where scalar perturbations are absent. One may find a suitable 
gauge - for instance the UFG, which is valid only during slow-rollover -, in which the back-reaction does not grow in 
time (because of cancellations of the infrared pieces). However, such a gauge seems inconvenient since the inflaton 
cannot be used as a clock in the subsequent stage of coherent oscillations. We feel that it is more proper to ask what 
the time behavior of back-reaction is in a frame which can be regularly continued rather than the opposite. Indeed, 
one of the advantages of the UCG is the possibility of continuing the calculation through the oscillatory phase of the 
inflaton [4(j, which is not easy to do in the UFG. It would be interesting to check if the same conclusions are obtained 
for the non-local quantity used in [4l| . 
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